Anisotropic upper critical field and a possible Fulde-Ferrel-Larkin-Ovchinnikov 
state in a stoichiometric pnictide superconductor LiFeAs 
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Measurements of the temperature and angular dependencies of the upper critical field H C 2 of a 
stoichiometric single crystal LiFeAs in pulsed magnetic fields up to 50 T were performed using a 
tunnel diode resonator. Complete H^(T) and H&%T) functions with H\l (0) = 17± 1 T, H^ 2 C (0) = 
26 ± 1 T, and the anisotropy parameter jh(T) = H^f/H'^ decreasing from 2.5 at T c to 1.5 at 
T<T C were obtained. The results for both orientations are in excellent agreement with a theory 
of H C 2 for two-band pairing in the clean limit. We show that H^(T) is mostly limited by the 
orbital pairbreaking, whereas the shape of i? c2 c (T) indicates strong paramagnetic Pauli limiting and 
the inhomogeneous Fulde-Ferrel-Larkin-Ovchinnikov (FFLO)state below Tf ~ 5 K. 



There are only few stoichiometric iron-based com- 
pounds (Fe-SCs) exhibiting ambient-pressure supercon- 
ductivity without doping. Among those LiFeAs is 
unique because of its relatively high T Q = 18 K, [l| as 
compared to LaFePO (T c = 5.6 K) [2] and KFe 2 As 2 
(T c = 3 K) [3j. The absence of doping- induced dis- 
order leads to weak electron scattering, low resistiv- 
ity, p(T c ) « 10 /iSlcm Q and high resistivity ratio, 
RRR = p(300K) /p(T c ) > 30 0, Q. These parame- 
ters differ significantly from those of most Fe-SCs for 
which superconductivity is induced by doping, for ex- 
ample, Ba(Fei_ a; T :E _) 2 As 2 @,Q,(Bai_ a: K 2; )Fe 2 As 2 i and 
BaFe 2 (Asi_ a; P a: ) 2 [8J. With the highest T c among stoi- 
chiometric Fe-SCs, negative dT c /dP Q, tetragonal crys- 
tal structure [HQ and the absence of antiferromagnetism 
[Toj | , LiFeAs serves as a model of clean, nearly optimally- 
doped Fe-SC 0]. Because of very high H c2 of Fe-SCs, 
they may also exhibit exotic behavior caused by strong 
magnetic fields, for example, the Fulde-Ferrel-Larkin- 
Ovchinnikov (FFLO) state in which the Zeeman split- 
ting results in oscillations of the order parameter along 
the field direction [ll|. Thus, measurements of H c2 (T) in 
stoichiometric LiFeAs single crystals can reveal manifes- 
tations of pairing in the clean limit [12j for which the 
FFLO state would be least suppressed by doping-induced 
disorder 11| as compared to other optimally doped Fe- 
SCs. 

Measurements of the upper critical fields parallel 
and perpendicular (H^f ) to the crystallographic c— axis 
in many Fe-Sc have shown several common trends [6, 
SEMI- Close to T c where H c2 is limited by orbital 
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pairbreaking, the anisotropy para meter y # = H^f/H'^ 



(Pc/Pab) 1/2 above T c 
more isotropic 




As T decreases, H c2 (T) becomes 
consistent with multiband 



ranges between 1.5 and 5 in agreement 



with the anisotropy of the normal state resistivity — 



miring scenarios and the behavior of H c2 in dirty MgB 2 
il], yet opposite to clean s ++ MgB 2 single crystals [29j j . 
However, the more isotropic H c2 at low T can also re- 
sult from strong Pauli pairbreaking for H||a6 since the 
observed H c2 on many Fe-SCs significantly exceeds the 
BCSparamagnetic limit H p {T) = IMT C {K) [13, El III - 
I27I |3G |. Thus, measuring H c2 in LiFeAs can probe the 
interplay of orbital and Pauli pairbreaking in the clean 
pairing limit at high magnetic fields. These measure- 
ments are also interesting because magnetic fluctuations 
may contain significant ferromagnetic contribution which 
may lead to triplet pairing [3ll ]. Experimentally, vortex 
properties of LiFeAs were found to be very similar to the 
supposedly triplet Sr 2 Ru04 [32(, although NMR studies 
suggest singlet pairing [33j ]. Triplet superconductors can 
exhibit unusual response to magnetic field [34| . and, in- 
deed, candidate materials show pronounced anomalies, 
as observed in UPt3 H 36| and Sr 2 Ru04 [33] • Surpris- 
ingly, our measurements show that normalized H^f of 
LiFeAs matches quite closely that of Sr 2 Ru04. 

We present the measurements of the complete H — T 
phase diagram of LiFeAs in pulsed magnetic fields up to 
50 T, and down to 0.6 K using a tunnel diode resonator 
(TDR) technique. We found that H^ 2 C {T) shows rapid 
saturation at low temperatures, consistent with strong 
Pauli pairbreaking. Similar conclusion was reached from 
torque measurements [4(j. Our data can be described 
well by a theory of H c2 for the multiband s pairing 
in the clean limit 38], which also suggests the FFLO 
state in LiFeAs for iJ_l_c below 5 K. Previous measure- 
ments of H c2 in LiFeAs were performed at relatively low 
fields [1 [3^ , thus not allowing to reveal the spin-limited 
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FIG. 1. (Color online) TDR frequency change for increasing 
pulsed magnetic field, appied in two orientations, H || c and 
HA-C, shown for two temperatures for sample A. The defini- 
tion of H C 2 is shown as the intersection of two straight lines 
below and above the transition. 



behavior at low T. The only reported high-field mea- 
surements associate -ff c2 with the disappearance of irre- 
versibility in torque measurements Ref. |40j . The authors 
supported this association by comparing with the specific 
heat data. However, in our opinion, the irreversibility 
field may underestimate the true H C 2(T) and have differ- 
ent temperature dependence due to depinning of vortices. 
It may also have significant (cusp like) angular variation, 
which would be particularly important for torque mea- 
surements that rely on the finite angle between magnetic 
moment and field. Related complications were discussed 
in high— Tc cuprates [ill ]. 

Single crystals of LiFeAs were grown in a sealed tung- 
sten crucible using Bridgeman method and placed in am- 
poules. Immediately after opening, samples were cov- 
ered with Apiezon N grease, which provides some degree 
of short-term protection The samples were cleaved 
and cut inside the grease layer to minimize exposure to 
the air. The two studied samples had dimensions of 
0.6 x 0.5 x 0.1 mm 3 (sample A) and 0.9 x 0.8 x 0.2 mm 3 
(sample B). Superconducting transition temperature for 
both samples was T c = 17.6 ± 0.1 K (more than 10% 
higher than T c = 15.5 K of Ref. lioh . Dynamic mag- 
netic susceptibility x wa s measured with 190 MHz (sam- 
ple A) and 16 MHz (sample B) TDR jH). The magnetic 
field was generated by a 50 T pulsed magnet with a 11 
ms rise time at Clark University. A single-axis rotator 
with a 0.5° angular resolution was used to accurately 
align the sample with respect to the c— axis (see inset in 
Fig. [2ja)). The data have been taken for each orienta- 
tion at temperatures down to 0.66 K. The normal state 
data at 25 K have also been taken for both orientations 
and subtracted. Measured shift of the resonant frequency 



A/ oc x 1421; thus exhibits a kink at ff c2 where London 
penetration depth diverges and is replaced by the nor- 
mal - state skin depth. Thus, barring uncertainty due to 
fluctuations, it is probing a "true" upper critical field. 

Fig. [T] shows the change of the resonant frequency as 
a function of H for sample A for two field orientations 
and two temperatures. From many such traces, both 
H^f and H^ were determined as shown in Fig. Q] 
and are plotted in Fig. [2] Figure EJa) compares our 
H C 2 data on samples A and B with the previous trans- 
port 0, |H, 43 1 and torque measurements [ipj . Fig- 
ure [2ja) also shows the behavior expected from the or- 
bital Werthamer-Helfand-Hohenberg (WHH) theory [IH 
with fforb(O) = 0.69T c \dH c2 /dT\ Tc , the single-gap BCS 
paramagnetic limit, H^ cs = 1.84T C = 32.2 T, as well 



as Hi 



= 34.7 T and H^ 2 = 20.4 T calculated with 



Ai(0)/T c w 1.885 and Ai(0]/T c « 1.111 reported for 
the same samples in Ref. 4]. Clearly, the observed 
H C 2(T) exhibits much stronger flattening at low temper- 
ature compared to the orbital WHH theory. Inset in 
Fig. [3£a) shows the dependence of H c2 on the angle be- 
tween H and the ab plane at 0.66 K where H^f is defined 
at a maximum of H C 2((p) = -ffj 2 + {H^f — H^ ) cos tp de- 
picted by the solid line. 

We analyze our H C 2 (T) data using a two-band theory, 
which takes into account both orbital and paramagnetic 
pairbrcaking in the clean limit, and the possibility of the 
FFLO with the wave vector Q(T, H ) . In this case the 
equation for H C 2 is given by [38| . 
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Here Q(T,H) is determined by the condition that 
H C 2(T,Q) is maximum, a\ = (Ao + A_)/2w, a 2 = 
(A - X-)/2w, A_ = A n - A 22 , A = (X 2 _ + 4A 12 A 2 i) 1 / 2 , 
w = AnA 22 — Ai 2 A 2 i, t — T/T c , and G 2 is obtained by 
replacing \fb — > y/rjb and q — > q^fs in G\ , where 
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h 2 v 2 



= e 2 /ex. (4) 



Here vi is the in-plane Fermi velocity in band I = 1,2, 
€4 = mf b /mf is the mass anisotropy ratio, cf>o is the flux 
quantum, /x is the magnetic moment of a quasiparticlc, 
An and A 22 are the intraband pairing constants, and Ai 2 
and A 2 i are the interband pairing constants, and a w 
0.56cva/ where the Maki parameter qm = H°J^ / 'V2H p 
quantifies the strength of the Zeeman pairbreaking. The 
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FIG. 2. (Color online) (a) H o2 (T) for H±c (solid symbols) 
and // || c (open symbols). Blue circles and red squares cor- 
respond to samples A and B, respectively. For comparison 
we show the literature data determined from the resistiv- 
ity measurements with mid-point criterion: ( mag enta) tri- 
angles [U, (green) rhombi [3!|, (brown) stars [4a |. Torque 
data are shown by (grey) pentagons Dashed lines is the 
WHH H c2 {T). Inset in (a) shows H c2 (<p) at 0.66 K where 
the solid line is H c2 (ifi) = H| 2 C + (H^ - H^)cosip. (b) Fit 
of the experimental data to H C 2(T), Q(T) and 7h(T) (solid 
lines) calculated from Eq. {1} for the parameters given in the 
text. The FFLO wave vector Q(T) is plotted in the units of 
40nkBT c gi/hvi, and the inset shows the anisotropy parame- 
ter 7ff(T). 



factors gi = l + An + |Ai 2 and g 2 = l + A 22 + |A 2 i| describe 
the strong coupling Eliashberg corrections. For the sake 
of simplicity, we consider here the case of ei = e 2 = e for 
which H^f is defined by Eqs. ([1]) and ^ with g\ — g 2 
and rescaled q —> qe -3 / 4 , a — > ae _1//2 and Vb — > e 1 / 4 -^ 
in Gt and -> e 1 / 4 ^ in G 2 [H. 

Figure [U[b) shows the fit of the measured ff c2 (T) to 



0.25, 
tion 



i] = 0.3, a 
(P) describes 



0.35, and e : 



0.128. Equa- 
and jh(T) = 



b||(T)/Ve6_L(T) where 6||(T) and 6 X (T) are the solu- 
tions of Eq. ([1]) for H\\c and if-Lc, very well. The 
fit parameters are also in good quantitative agreement 
with experiment. For instance, the Fermi velocity v\ — 
(gik B T c /h)[8TT(t)ob±(0)/Hl2{0)} 1/2 can be expressed from 
Eq. ^ in terms of materials parameters and b±(0) = 
0.314 calculated from Eq. Q. For T c = 17.8 K, #J 2 C (0) = 
18.4Tand.g = 1.5 for Ai 2 = 0.5, we obtain vi = 1.12xl0 7 
cm/s, consistent with the ARPES results [10| . 

Several important conclusions follow from the results 
shown in Fig. HJb). First, contrary to the single-band 



Ginzburg-Landau scaling, 7^ 
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the anisotropy 



Eq. ([1} for ,s ± pairing with An = A 22 = 0, Ai 2 A 2 



parameter 7h(T) decreases as T decreases. Not only is 
this behavior indicative of multiband pairing 281 ] . but it 
also reflects the significant role of the Zeeman pairbreak- 
ing in LiFeAs given that ay = a/y/e — 0.98 for H±c 
is close to the single-band FFLO instability threshold, 
awl [38]. In this case 70 (T) near T c is determined by 
the orbital pairbreaking and the mass anisotropy e, but 
as T decreases, the contribution of the isotropic param- 
agnetic pairbreaking increases, resulting in the decrease 
of 7h (T) . Another intriguing result is that the solution of 
Eq. ([T|) shows no FFLO instability for H\\c, but predicts 
the FFLO transition at T < T F sa 5 K for H\\ab. The 
FFLO wave vector Q(T) = ^k B T c q{T)b 1 / 2 {T)g 1 /hv 1 
appears spontaneously at T = Tp « 5 K where the 
FFLO period I = 2-k/Q = hv 1 /2k B T c g 1 q(T)b 1 / 2 {T) di- 
verges and then decreases as T decreases, reaching £(0) = 
7rV3i<?(0)& 1/2 (0) ~ 9£ at T = 0. Here q(0) = 0.656, 
6(0) = 0.126, and £ = hv 1 /2nk B T c ~ 7.3 nm, giv- 
ing £(0) ~ 65.6 nm for the parameters used above. 
The period £(0) is much smaller than the mean free 
path, £ m f p ~ 550 nm, estimated from the Drude for- 
mula for an ellipsoidal Fermi surface with e = 0.128, 
vp = 112 km/s, m a b equal to the free electron mass, and 
p(T c ) — 10/iflcm. Notice that p{T c ) may contain a signif- 
icant contribution from inelastic scattering, so the mean 
free path for elastic impurity scattering which destroys 
the FFLO state [H| is even larger than l m f p . Therefore, 
the FFLO state found in our calculations may be a realis- 
tic possibility verifiable by specific heat, magnetic torque 
and thermal conductivity measurements. 

Finally, we compare LiFeAs with other supercon- 
ductors, especially those for which i? c2 is clearly lim- 
ited by either orbital or Zeeman pairbreaking mecha- 
nisms. Shown in Fig. [3] are the plots of the normal- 
ized H c2 (T)/T c H' c2 as functions of T/T c for the H||a6 
orientation where the Zeeman pairbreaking is most pro- 
nounced. Here H' c2 = \dH C 2/dT\T-tT c an d our data are 
shown by the thick solid black line, whereas the liter- 
ature data are shown by symbols. The reference ma- 
terials include i/ c2 for: LiFeAs Pauli- limited [45] ] 
organic superconductor K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br 



4 



0.6 



rr 0.4 



X 



N 
c 

X 



0.2 - 



0.0 



aJ 1 


1 I 


1 1 1 


LiFeAs ^Njl 
















-A-Fe, „Te„ Se„, 

111 0.6 0.4 




" — *— Ba K FeAs 

0.6 4 2 




BaFe ,. B 4 C °0, 6 AS 2 


LiFeAs 




-<- NdFeAsO, ^ 
■ -B-KFe 2 As 2 


this work 

Kurita et al. 




-»-K-(ET) 2 Cu(N(CN) 2 )Br 


-y- Sr 2 Ru0 4 (H||a) 


-©-CeColn 5 

i , i 


-O-NbTi aR 

■ 



0.0 



0.2 



0.4 



T/T 



0.6 



0.8 1.0 



FIG. 3. (Color online) (H c2 (T)/T c )/\dH c2 /dT\ Tc vs. T/T c in 
the H-Lc orientation. Black solid line is our data in compari- 
son with several Fe-SCs as well as other exotic superconduc- 
tors and conventional NbTi, all shown in the legend. 



[46| ; heavy fermion CeCoIn 5 14 71 : optimally-doped iron 
pnictides, Ba(Fei_ 2; Co a ;)2As2 [18[ and Ba x Ki_ x FeAs2 



[20j as well as iron chalcogenide Fe(Se,Te) [27|. Con- 
ventional NbTi is also shown by open pentagons jijj . 
Remarkably, scaled data obtained on crystals with dif- 
ferent T c s and by different measurements (this work and 
Ref. [loj]) are very similar indicating intrinsic behavior 
of LiFeAs, namely, that it is indeed closer to the para- 
magnetic limit. Notably, the data for LiFeAs lay below 
other Fe-SCs, except for the highest purity (RRR ~ 87) 
KFe2As2 [49j . On the other hand, our data appear above 
CeCoIn 5 , believed to be mostly Pauli limited [47J. Inter- 
estingly, the data for LiFeAs stay almost on top of the 
H C 2(T) for Sr2Ru04, in which limiting of H C 2 proceeds 
in a very unusual manner, leading to the formation of 
the second superconducting phase [37j. Given that vor- 
tex dynamics in these two materials is also similar [32| . 
the coincidence of the H C 2(T)/T c H' c2 curves is worth of 
further exploration. 

Summarizing, full - temperature range experimental 
(T) and H^f(T) deviate significantly from the single- 
band WHH behavior but are in excellent agreement 
with the theory of H C 2 for the pairing in the clean 
limit. Our results indicate Pauli-limited behavior and 
the FFLO state below 5 K for H±c. 
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